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Abstract. This paper studies the Casimir effect due to fractional mass- 
less Klein-Gordon field confined to parallel plates. A new kind of bound- 
ary condition called fractional Neumann condition which involves van- 
ishing fractional derivatives of the field is introduced. The fractional 
Neumann condition allows the interpolation of Dirichlet and Neumann 
conditions imposed on the two plates. There exists a transition value 
in the difference between the orders of the fractional Neumann condi- 
tions for which the Casimir force changes from attractive to repulsive. 
Low and high temperature limits of Casimir energy and pressure are ob- 
tained. For sufficiently high temperature, these quantities are dominated 
by terms independent of the boundary conditions. Finally, validity of 
the temperature inversion symmetry for various boundary conditions is 
discussed. 

PACS numbers: ll.10.Wx 



1. Introduction 

Applications of fractional calculus, in particular fractional differential 
equations, in transport phenomena in complex and disordered media have 
attracted considerable attention during the past two decades [H [2J El HI El E] ■ 
However, the use of fractional calculus in quantum theory is still very lim- 
ited. Recently, generalization of quantum mechanics based on fractional 
Schrodinger equation has been considered by several authors [7] [51 I§1 flu] [TT] . 
In quantum field theory, fractional Klein-Gordon equation |12 [ ll3 [ rH lll5[ ll6] 
and fractional Dirac equation [171 [18] were introduced several years ago, 
but further studies on these topics are scarce. It was only lately that the 



Key words and phrases. Casimir energy, fractional Klein-Gordon field, fractional Neu- 
mann conditions, temperature inversion symmetry. 

1 



2 



C.H. EAB\ S.C. LIM 2 , AND L.P. TEO 3 



canonical and stochastic quantization of fractional Klein-Gordon field and 
fractional Maxwell field have been carried out [191 [20} \2l\ [22] . 

In this paper, we shall consider another aspect of fractional Klein-Gordon 
field, namely, the Casimir effect associated with such a field. This work 
is partly motivated by the recent advances in cosmology, in particular the 
solid evidences for accelerated expansion of the universe [23l [24], [25J [26] . 
which have rekindled considerable interest in Casimir effect [27, 28J. Casimir 
energy in extra space-time dimensions [29] has been proposed as a possible 
candidates of dark energy [30] that is responsible for the accelerated cosmic 
expansion. However, in this paper, we shall not deal directly on the link 
between Casimir energy and the dark energy. Instead, we shall study the link 
between the possible repulsive nature of the Casimir force and the general 
boundary conditions associated with fractional massless Klein-Gordon field. 

In most consideration of Casimir energy between a pair of parallel plates, 
the boundary conditions employed are either of Dirichlet type or Neumann 
type for both of the plates. A less common pair of parallel plates has been 
suggested by Boyer [31], with one of them perfectly conducting and the 
other infinitely permeable. Boyer was able to show in the context of ran- 
dom electrodynamics that for such a set-up the resulting Casimir force is 
repulsive. It is possible to show that this unusual pair of plates necessi- 
tates mixed boundary conditions, with the Dirichlet condition for the per- 
fectly conducting plate, and Neumann condition for the infinitely permeable 
plate. Recently, this result has been derived by several authors using the 
zeta function method for scalar massless field at zero temperature [32] and 
finite temperature [33] , 

Since this paper studies Casimir effect associated with fractional Klein- 
Gordon field, it is not unnatural for one to consider the fractional generaliza- 
tion of Neumann conditions involving fractional derivatives. We shall study 
how repulsive Casimir force due to the fractional massless Klein-Gordon 
field can arise under a new type of boundary conditions, the fractional de- 
rivative boundary conditions (or fractional Neumann conditions). We show 
that such conditions allow interpolation between the ordinary Dirichlet and 
Neumann conditions. 

This paper is organized as follows. In next section we first recall some 
basic facts about fractional Klein-Gordon field at zero and finite tempera- 
ture. In Section 3, we derive the partition function and free energy between 
parallel plates associated with the fractional massless scalar field at posi- 
tive temperature using the generalized thermal zeta function regularization 
technique. We show that the Casimir force associated with the massless 
fractional scalar field can change from attraction to repulsion as the order of 
the fractional Neumann conditions imposed on the parallel plates is varied. 
Finally we obtain the low and high temperature limits of various physi- 
cal quantities such as free energy and pressure. The temperature inversion 
symmetry [33j [Ml E3 HH EH [38l E2] will also be discussed. 
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2. Fractional Klein Gordon Field 

In this section, we recall briefly some basic theory of fractional derivative 
fields. Let us consider the Euclidean scalar field </>(x, t), x£ M> D , t £ 1R with 
the following Lagrangian 

(2.1) L = ^(x, t )A(-A)0(x,i), 

where A = d 2 + Ylf=i ®] ^ s the (D + l)-dimensional Euclidean Laplacian 
operator, and A(— A) is a pseudo-differential operator [40J. In order to con- 
sider A(— A) of fractional order which contains the fractional powers (— A) Q 
of Laplacian operator, we need to define the Riesz fractional derivative and 
integral [UJ in order to give these operators a precise meaning. For a test 
function in Schwartz space (or a tempered distribution) g, the Fourier trans- 
form of — (Ag)(x) satisfies — A<?(£) = \£,\ 2 g(£,)- This can be generalized to 
fractional power of Laplacian operator. For our purpose, it is sufficient to 
consider only the real fractional powers. For a G M \ {0}, and Schwartz 
functions g we define 



(2.2) (-A)-°/ 2 5 (x) = (ier/(e)) v (x) 



\l a g(x), a>0, 
VD~ a g(x), a < 0. 

The operators I Q and D Q defined in (|2.2p for a > are called respectively the 
Riesz fractional integral operator and Riesz fractional differential operator. 
We have B a I a g = g and I a Pg = l a+l3 g , a > 0,(3 > for "sufficiently 
good" functions g. 

A(— A) in (|2.ip can be expanded in a power series ^ ■ Cj(— A) J , and it can 
be regarded as a differential operator of infinite order of derivatives. From 
the Lagrangian field theory with higher order derivatives [20, 42J one gets 

3 3 

and by summing up the series gives the nonlocal field equation A(— A)</>(x, t) = 
0. Nonlocal field theory with A(— A) = (—A + m 2 ) a , a > as the fractional 
Klein-Gordon operator has been considered by several authors |12l [T3| HU 
HS| [T6| [T9l [20l [21] [22] . Higher derivative field theories involving propagator 
of the form (k 2 + m 2 ) _n , n > 1 were first used by Pais and Uhlenbeck [33] to 
obtain a regularized theory without ultraviolet behavior. Fields with such 
propagators result in either theories with ghost states that require a Hilbert 
space with indefinite metric, or nonlocal theories without ghost states. 

Here we give some remarks on the motivations for introducing fractional 
derivative fields. Field theories with nonlocal Lagrangian of the type (|2.ip 
with nolocality due to kinetic terms have attracted considerable interest. 
For examples, nonlocal kinetic term plays an important role in the (2+1)- 
dimensional bosonization [441 145] : it also arises in effective field theories 
when some degrees of freedom are integrated out in the underlying local 
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field theory [461 l4"Tj . One also expects fractional derivative quantum fields 
to play an important role in quantum theories of mesoscopic systems and 
soft condensed matter which exhibit fractal character. Such argument can 
be extended to quantum field theories in fractal space-time [48 \ [49] . 

Canonical quantization of nonlocal scalar fractional Klein-Gordon field 
has been considered by Amaral and Marino [19], and Barci, Oxman and 
Rocca [20] . Free relativistic wave equations with fractional powers of D'Alembertian 
operator were studied by several authors [HI [TU [I5j [16]. Stochastic quan- 
tization of fractional Klein-Gordon and fractional abelian gauge field has 
been considered by Lim and Muniandy [21], and finite temperature frac- 
tional Klein-Gordon field is considered in a recent work |22j . The two-point 
Schwinger function of the Euclidean fractional Klein-Gordon field is given 



For the Euclidean fractional Klein-Gordon field at finite temperature T = 
1/(3 satisfying the periodic condition cp(x,y, z,0) = <f)(x,y, z, /3) , the two- 
point Schwinger function becomes 



where k n = (k, u n ), uj n = 2nir/ (3 and = k 2 +w^. The two point Schwinger 
functions for the massless field are given by (|2.3p and (|2.4p by putting m = 0. 

In the next section, we shall carry out the computation of Casimir en- 
ergy associated with the massless fractional Klein-Gordon field confined 
between two parallel plates imposed with fractional Neumann boundary 
conditions. The thermal zeta function technique will be employed in our 
calculation. Zeta function method was introduced as regularization pro- 
cedure in quantum field theory about two decades ago [501 ED E2]- Basi- 
cally the zeta function technique involves three steps. In the case for scalar 
massless fractional Klein-Gordon field they are: (I) Determination of the 
eigenvalues A of (— A) a with appropriate boundary conditions, hence the 
spectral zeta function C(-A) a ( s ) = 0--0 Analytic continuation of 

the zeta function C(-A) Q ( S ) to a meromorphic function of the entire com- 
plex plane. (Ill) Evaluation of det(— A) a in terms of C(-A) a ( s )> that is, 



det(— A) a = exp C(_a)«(^)J ■ ^ or simplicity, the computation will be car- 
ried out for scalar massless fractional Klein-Gordon field. However, one can 
mimic the electromagnetic field by the scalar massless field with the two 
transverse polarization states of the former taken care of by multiplying the 
end results by a factor of two plus some minor modifications on the possible 
eigenmodes of the field. In this way, we can compare our results to some 
other established results. 



by 



(2.3) 




(2.4) 
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3. Free Energy of Massless Fractional Klein Gordon Field at 

Finite Temperature 

We first assume that the fractional Klein-Gordon field 0(x, t) is inside 
a .D-dimensional space f2 which is a rectangular box f2 = [0, Li] x . . . x 
[0,L£i_i] x [0,d] such that d <C 1 < i < D — 1. At the end, we 
let Lj, 1 < i < D — 1 approach infinity to obtain space between the two 
hyperplanes xjj = and xd = d in We want to consider massless 

fractional Klein-Gordon field confined in the region Q, and maintained in 
thermal equilibrium at temperature T = 1/(3. As usual [53], we impose 
periodic boundary condition with period (3 on the imaginary time, i.e. 

0(x,i) = 0(x,i + /3), VteR. 

The Helmholtz free energy of the system is then given by the equation 

1 







logZ, 



where Z is the partition function defined by 

(3.1) Z = J V[4>] exp (-± J J>(x, t)*(-A) a <t>(x, t)\ d D xdt. 

BC 

Here BC denotes boundary conditions on the field <fi. We impose periodic 
boundary conditions in the directions of x±,... ,xd-i- In the direction xd, 
we can consider different boundary conditions, among them are the Dirichlet 
boundary condition with 



(x,0,t) = ^(x,d,i) = 0, 



V x € 



which corresponds to perfectly conducting plates in the case of electromag- 
netic field; the Neumann boundary condition with 







dxD 



(x, x D ,t) 



d 



x D =0 



dxo 



(x,x D ,t) 



o.. 



V x G 



D-D-l 



t € 



X£)=d 



which corresponds to infinitely permeable plates in the case of electromag- 
netic field; and the mixed boundary condition with 

d 



(x,0,t) = 0, 



dx 



D 



(x,XD,i, 



o. 



V x e 



d-D-1 



X£)=d 



which corresponds to Boyer's setup (namely one plate is perfectly conduct- 
ing while the other infinitely permeable) in the case of electromagnetic field. 
Since we consider the Casimir effect associated with fractional massless 
Klein-Gordon field, one can consider the most general boundary conditions, 
namely the fractional boundary conditions 



(3.2) 



dx^ 



(x,x D ,t) 



0. 



x D =0 



9x d' 



{k,x D ,t) 



0.. 



XD- 
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where x-,^ £ [0)1]- Here we use the definition of fractional derivative in 
terms of Fourier transform: 

where 

(±ifc) Q = |fc| Q e ±1 ^sgn(A:) 

and f(k) is the Fourier transform of /. 

When x = A 1 = 0, one gets the Dirichlet condition for both the plates. On 
the other hand, when % = = 1, the boundary conditions for both plates are 
that of Neumann type. In the case with either x = 0, \i = 1 or x = 1> A 4 = 0> 
we have the Boyer type boundary condition. For values of (%, /i) other than 
the above values, we have fractional Neumann boundary condition for both 
plates. One can naively regard such boundary conditions as correspond to 
plates which are not perfectly conducting or infinitely permeable. 

Now we want to analyze the condition (|3.2p . If 

rp k {z) = Ae ikz + Be~ ikz 

are eigen-modes on the z = xd direction, the requirement (|3.2|) is equivalent 
to 

A(ik) x +B(-ik) x = 0, 
A{ikYe ikd +B{-ikYe~ ikd = 0. 
From these equations, we find that 

(3.3) B = -e inx A, 

(3.4) 2iAk x+ v e -^ sin (kd - |(x - m)) = 0. 
From (|3.4p . we find that the value of k has to be 

Together with (|3.3p . the eigen-modes in z direction are given by 

where rj = /i — x- We have the following specific cases: 
x = n = o, ip = 0, ip n = -tp- n ; 
x = fj, = i s ^ n = ijj_ n - 

X = 0,H=1, lp n = -Tpl-n, 
X = 1,/X = 0, Vn = ^-l-n; 

For all other cases, ip n ,n £Z are linearly independent. 
We let 5 X , M = N if ( X , //) = (0, 0) or (0, 1), S X , M = N U {0} if ( X , /x) = (1, 0) 
or (1,1) and S Xjfl = Z for all other cases so that {ip n (z) : n 6 S x ^} is 
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a complete set of linearly independent eigen-modes satisfying the condition 
(|3.2j) . Now it follows that the eigen-modes of the field </>(x, i) are 

[2mkixi\ (2-Kik D ^ix D -i\ [2mmt\ 
0k,n,m(x,z,t) = exp I — — — I . . . exp I ^— — - I Vw)exp I — - — I 

with k = (ki, ... , A;d-i) 6 Z d_1 , m G Z, n 6 <S XjA1 . 

As is well-known, up to a normalization constant, the path integral (|3.ip 
is equal to 

(3.5) Z m =l 11 11 J] 'A k , n , m j = [det(-Ar 1- ' " 



D-l /n i \2 



2irkj \ , fx ( 7] 




where 



The prime ' in ()3.5[) indicates the omission of X^ n) m = terms. We com- 
pute (|3,5p using zeta regularization [5H [551 ES]. Namely, we define the 
spectral zeta function 

(3-6) c«*» = E E E ' A k,n, m - 

Then 

logZ Q;XiM = ^C ;XlM (0). 

Obviously, Ca;x,Ai( s ) = Ca;M,x( s )- ^ n terms of the spectral zeta function, the 
Helmholtz free energy can be expressed as 

(3-7) F a . x>li = -^C, X ,M 

We are interested in the limit Li — > oo for all 1 < i < D — 1. In that case, 
instead of the free energy, we consider the free energy density 

(3-8) / a;XiM = ^^, where A = L 1 ...L D _ 1 . 

As usual, the pressure is related to the free energy by the thermodynamic 
formula 



(3.9) P c 



8V ) T \ 8(1 j rp 



In order to compute the spectral zeta function Co;x,At( s ); recall that the 
Epstein Zeta function is defined by (see e.g. [53] ) 

^ g27rin.h 

(3.10) Z B (s;ai,a 2 ;g,h) = V'— — ■ r=- ■ 

^ (ai(ni +5i) 2 + 02(^2 + 52) 2 ) s 
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Here as usual, the ' over the summation means that the term n = is 
omitted when g = 0. The function Z E (s; 0,1,0,2; g, h) satisfies the functional 
equation (see e.g. [54], page 6) 

(3.11) 

-27rig.h / 11 \ 

7r- s T(s)Z E ( S ; a x , a 2 ; g, h) = - ^"^(l - s)Z E 1 - s; -, -; h, -g 

y/aia 2 V ai a 2 J 

In the following, we carry out the computation of Ca X ^(0) f° r various 
boundary conditions. 

3.1. The case X ^ V and ( X , n) + (0, 1), (1, 0). 
This corresponds to the boundary condition 



QX 

—<f>(x.,x D ,t) 



dx D 



°' TTjr4>(x,x D ,t) 

00c J 



= 0, 

X£)=d 



x D =0 D 

< x < 1, < /j, < 1. 

In this case, 77 7^ 0, ±1. The zeta function Cq;x,m is given explicitly by 

c™w-EE E (E(^) 2 + e("-f)) 2 + (^) 

neZmeZkez 15 - 1 \j=i v J 7 v ' 

To simplify notation, let 

/7r\ 2 f 2ir\ 2 77 

ai= y ' ° 2= ' c= 2- 

As Lj — » 00 for all 1 < i < D — 1, 
(3.12) 



(2 7 r) D - 1 r , f^_ 2 7r [k 2 + ai(n- c) 2 + a 2 m 

27r (D-l)/2^ ,00 



( 27 r)^r(5-i) 



fc D - 2 (jt 2 + iy as dkx 
1 



E 



, (ai(n-c) 2 + a 2 m 2 ) a — K^- 1 )/ 2 ) 

A r(as-^i) / (X>-1) ,\ 
- -Z E \as ;ai,a 2 ;g,h , 



(4tt)( d - 1 )/ 2 r(as) 
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with g = (— c, 0), h = 0. From the functional equation (|3.11|) . we find that 
(3.13) 

A7T 2as-D r(^±i_ as )^ fD + 1 11 



Ca;xA S ) = Ta \(D 1V2 / TV \ Ze O aS > ~ ' _ ' h ' _ § 

(47r) ( .- u_i i/ i A /a 1 a 2 T(as) \ 2 a± a 2 

A7T 2as-D T (D±l _ as \ e 2™ c 

E 2J-, 



WP-W^/aiaZ T(as) ([„2/ ai ] + [ m 2 /a2])[ p + i)/2]- 

This gives 

C' (o) = 1 2 j V V- 

Sa^./xV"^ oD-l„(3D-l)/2 /77777T Z_> , 



2 D-l 7r (3D-l)/2^ i Z_. ^ ([n 2 /ai] + [m 2 /a2]) (D + l)/2 



r (^) 



2ar i)/2 ^P + l) 



2 D-l vr (3D-l)/2^ T ^ 

, 2 ( D +!)/2 cos(27rnc) , cos(27rnc) 

n D + l + A,Z l([n 2 /ai] + [m 2 /a2])(D+ l)/2^ 

where Cr( s ) is the Riemann zeta function. Define 

£ — — J_ 

7T/3 27rY Ol 



In terms of £, the free energy density f|3.8j) is equal to 
(3.14) 

2ar (jg) vr^ / 1 « cos^m?) 

/ OT ~p U UD + 1) + ^ 



71=1 



2e D+1 E E 



cos(7rnr}) 



(™ 2 + (2<n) 2 )( D + 1 )/ 2 y ' 
and the pressure (|3.9|) is given by 
(3.15) 



D + l 



_ 2ar(^)vr— / 0+ p ^ cosfrm,) 



71=1 



n 1 



2 .D+i cos(7rn>7)(m 2 - £>(2<n) 2 ) \ 

^i^i ( m2 + (2vr£n) 2 )( D + 3 )/ 2 J' 



By using the formula 9.622 in |57] . 

(-l) n - 1 2(2n)! ^ cos (27rA:a;) 

where i?2n is the Bernoulli polynomial of order 2n. In particular, when D = 
3, using B A (x) = x 4 - 2x 3 + x 2 - 1/(30) and Ch(4) = -vr 4 J B 4 (0)/3 = vr 4 /90, 
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we find that the free energy density and the pressure are given respectively 
by 

/ fi«-4 2 oooo / \ \ 

a I 7r u 4 7r /rj\ 2t 4 cos(7rn?7j \ 



a;x ' M d 3 I 45 24 4 V2,/ s ^ ^ (m 2 + (2<n) 2 ) 2 

\ m=l n=l v v ' ' 



P - I Z__e! 4. !l! r (T[m a^m V V cos^nr?)^ 2 -3(2<n) 2 ) 
" d 4 1 45 8 v 2 / ^ ^ (m 2 + (2<n) 2 ) 3 

\ m=l n=l v ^ ' ' , 



3.2. The case x = M = [Dirichlet Boundary Condition] 

In this case, the zeta function C«;0,o (|3.6p is given by 

(3.16) 



EE E E(^) 2 + (t) 2+ (t) : 



Ca;0,o(s) = 

n=lmgZkeZ 13 - 1 \i 



Using the same method as in Section 13.11 we obtain 
(3.17) 

c ( S )= A r ( ag -^) yy 1 

W;0,0W (47r)(D -l)/2 r(as ) ( 2 + 2)a S -[(D-l)/2] 

v ' n=l mgZ v ' 

^rfas-Ari) / / D-l 

as — ; ai, a2 



2(4vr)( D - 1 )/ 2 r(as) 



- 2a [ i D - 1)/2] - as ( R (2as - (D - 1))J 
^ /^ 2 « S - D r (^±i - as) /D + l 11 



— „ as; 



2(4vr)( D - 1 )/ 2 r(as) I y/aiaj \ 2 ' a x ' a 2 



Here we have used the functional equation 



r © C/?(s) = rfr ( V) Cr{i - s) 

for Riemann zeta function. The first term in (|3.17j) is half of the term (I3.13j) 
with rj = 0. Therefore, we find that the free energy density and the pressure 
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are given respectively by 
(3.18) 

/a;0,0- ^ ^ Cfl(^ + 1)+ (2yr)D+1 - 2 ^ r( D±l^ CrW 

oo oo ^ \ 



(™ 2 + (2<n) 2 )( D + 1 )/ 2 j ' 



^a;0,0 TnT^ K ChC^ + 1 



d D+i y > (2tt) d + 1 

cos(7rn?7)(m 2 — D{2ir^n) 2 ) 



!i ^ (™ 2 + (2vr£n) 2 )( D + 3 )/ 2 ^ ' 
In particular, when Z) = 3, 

/«;0,o = -^3 I -45- + too + 4w ^ 2^ 2^ (m2 + (2^)2)2 - — Cr(3) I , 

\ m=l n=l v v 7 ' / 

and 

= ^Ar^_^ , ^ ^ m 2 -3(2^n) 2 \ 

a;0 '° 2d 4 I 45 240 ? ^ ^ (m 2 + (27r£n) 2 ) 3 i ' 

\ m=l n=l v v ' ' J 



3.3. The case x = ^ / 0, 1. 

This corresponds to the boundary condition 



(x,XD,t) 



x D =0 



(x,XD,t) 



0. 



0<X<1. 



xn=d 



In this case, rj = and the associated zeta function (|3.6p becomes 
(3.19) 

2 /2?rm x 2 



Ca;x,x( s ) 



E E 

(k,n,m)eZ D + 1 \{0} \j=l 



/ 27rkj\ 2 /irn 



which can be written as the sum of two terms 



E E (E(^) 

keZ 13 - 1 (n,m)GZ 2 \{(0,0)} \i=l V J 7 



2tt/v, \ ~ ^Trn\ 2 ( 2nm 



\ d 



+ 



+ E E 

keZ D-l\{ } V i=i 



/ 2irkj x 



V /3 
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The first term Ca-xx( s ) can De com P u ted as in Section I37L1 and the result is 
the same as ()3. 13|) with 77 = 0. For the second term Ca;y,v( s )' we wan t to 
verify in the following that it does not contribute to the free energy density. 
We have 

CM =2 E E i D f' 2 ' k ' 

(fci,...,fc D _ 2 )ez-D-2 fcD _ 1= i V j= i 



E E(^) I =««+«(•)■ 

(fei,...,fc £ ,_ 2 )6Z £, - 2 \{0} \i=i 
In the limit Lj — ► 00 for all 1 < i < D — 1, 




= 2 E EE 

(fcl,...,fc D _ 2 )eZ- D -2A: I 3_ 1 =l \j=l 



L 1 ! " ** 

/ 2irkj Si 2 



A/ 



. 2L i--- L ^-2 v / rf°- 2 k I Ikl 2 + f^Ezl V I 

«D— 1— -ImD— 2 v ' 



(2^-r(^) V V^-i 



(2tt)^ 2 r(Qs) ^ V Ld _ 



2vr(^ 2 )/ 2 L 1 ...L D _ 2 ( 2vr \ D - 2 ~ 2 « S r (as - 2=*-) 



( R (2as - (D - 2; 



(2vr)^- 2 V^-i/ r(aa) 

2QS _3(£^)±i D _ 2 _ 2as 

f 1 r — as )(r(D-1 



Therefore, 

3(D-2) + l , . n 9 

vr 2 L X ...L D _ 2 2vr \ 2 / D - 1 



v ><°>= — (i^rj ^_jc R (D-i) 

and the limit Xvaxi li ^, 00 {Y[{^) I (L\ . ..Ld_i)) vanishes. Similarly, the limit 
linii i ^ oo (y 2 / (0)/(Li . . . Ld-i)) = 0. Consequently, the contribution to the 
free energy density only comes from Ccrx,x'(^) ana " we ^ na - that the free 
energy density and the pressure in this case are given respectively by (|3.14p 
and (|3.15p by putting rj = 0. In particular, when D = 3, 

V + h + ^ T,E {m 2 + (2^)2)2 > 
m=l n=l v v 7 ' / 
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Q;x ' x d 4 I 45 240 ? ^ ^ (m 2 + (2vr£n) 2 ) 3 

V m=l n=l v , 



3.4. The case x = M = 1 [Neumann Boundary Condition]. 

In this case, the corresponding zeta function (|3.6p is given by 



, -ccs 

c,u = E E E'(e(^) 2 + (t) 2 + (T) 2 

n^mGZkeZ-D- 1 \i=l V J 7 ' 

It is easy to see that the sum of Co;i,i with Ca ; o,o f)3. 16j) gives Ca;x,x> X 7^ 0, 1 
(gllSl) - Therefore 

C;l J l(0)=C; x , X (0)-C;O,0(0)- 

We obtain from (l3TT8|) in Section O and ([3TT4"]) in Section O (with r? = 0) 
that in this case, the free energy density is given by 

^ U C«P + i)+ (27r)D+1 + 2V?r( D±l)^ C ^ 



2^ E E 



2 

1 \ 



( m 2 + ( 27r f„)2)P+l)/2 I ' 
m=l n=l v v ' ' / 

and the pressure P a -,i,i = P a -fl,o- hi particular, when D = 3, 



f a ( ^ 4 + ^ 4. 4.V V V 1 , ^\ \ 

f ^ l = ~2d? 720 ^ ^^KT(W ~ (fi( )| ' 

\ m=l n=l v v ' ' 



3.5. The case {xi A 4 ) = (1, 0) or (0,1) [Boyer Boundary Condition]. 

In this case, the corresponding zeta function (|3,6p becomes 

/ \ —as 

/o r \ 2 / / n \\2 /o \ 2 



c«iW=EE E s(^y + (5(-i))- 



+ 



n=imeZkez fl - 1 \i=i 
Observe that 



w W 4ee e (e(^) 2 ^H)H^) 2 
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Therefore, we can obtain the free energy density for this case by multiplying 
(|3.14p in Section [XT] by 1/2 and setting 77 = 1. This gives us 

U,o,i- jo K (r{d + 1)- ^ D+1 ( R (D + i) 

00 00 f i\n \ 

+ 2( D+1 V V ^ 

. . (™ 2 + (2vren) 2 )( D + 1 )/ 2 ) ' 

m=l n=l v v ' ' / 



and 



_ «r(^±i)^ / ^(i-2- g ) , 

fa;0,l- ^+1 K ChC-D + I) ( 27r )D+l CaC^ + 1) 

, 2 ,D+i V V (-l)"K-^(2<n) 2 ) \ 
^i^i (m 2 + (2<n) 2 )(^+ 3 )/ 2 )' 



When L> = 3, 



a ( 7T 6 £ 4 7 7T 2 oo oo (-1)™ 

(3.21) / a ;0,l = "^3 I ^ " g ^ + 4 ^ E E (m 2 + (2^)2)2 
\ m=l n=l x , 

_ a /A 4 , 7 ^ 2 , 2 4 ^ " (-ir(m 2 -3(2<n) 2 ) 
01011 2d 4 I 45 8 240d 4 ? ^ ^ (m 2 + (27r£n) 2 ) 3 

\ m=l n=l 



The results obtained for various boundary conditions can now be summa- 
rized in the following compact form: 

(3.22) 

g^oT 7T— /^ D+1 ^ . 1 ^ cos(7rnr?) 

/( w = dP ^ + 1} + (27pTT n D+i 

\ n=l 
I ^ ^ CQS(TO7?) r (f) t D A 

hh (- 2 + (2vrCn) 2 )(^)/ 2 + ^2^T(^f ^ ) ' 

_ f J X;/ ,ar(^)7r i ^ 1 m . g ^ cosing) 

^«;x.A» - d D+i I « + ^ (27r) c +! ^ n D + l 

V 71=1 

D+i A A cos(7rn>7)(m 2 - D(2<n) 2 ) \ 
^i^i (m 2 + (2vr^n) 2 )( D + 3 )/ 2 )' 

m=l n=l v v 7 ' / 
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where 



a 



uJ 



X,M 



1, 



0. 



if ( X ,fi) = (0,0), (0,1), (1,0), (1,1); 
else. 

if ( x ,m) = (1i1); 

if (x,A») = (0,0); 
else. 



Tracking back the derivation of formula (|3.22[) , we can also write the free 
energy density as (see (|3.13p ) 



(3.23) 



with g = (-r//2,0). 



a^adT{^) (D + l_ 1 1 _ 



2 a% ci2 



avrTr(^) n , 



2d D 



3.6. Casimir Energy of Electromagnetic field confined between par- 
allel walls. 

As is well known (see e.g. |53j), the Casimir energy of electromagnetic 
field in four dimensional space-time confined between two infinite parallel 
plates can be computed using almost the same setup as the massless scalar 
field with D = 3 and a = 1. More specifically, since there are two transverse 
polarization for electromagnetic fields, its free energy will be twice that of 
the massless scalar field. In the case when the two parallel plates are both 
perfectly conduction, except for the factor of two, it is almost equivalent to 
the Dirichlet boundary condition. However, as pointed out in [581 153| [59]. 
an additional 1/2 of the n = modes must be added. This amounts to the 
omission of the second term in (|3.17p . Therefore, the Casimir energy density 
for the fractional electromagnetic field is given by 

m=l n=l v v ' ' J 

in perfect agreement with the result obtained in [60] when a = 1. In Boyer's 
setup, where one plate is perfectly conduction and the other is infinitely 
permeable, the result for electromagnetic field should be twice the result for 
massless scalar field under Boyer's boundary condition. In fact, when D = 3 
and a = 1, twice of the formula (|3.2ip agree with the result obtained in [33] . 

By these comparisons with electromagnetic field, one can provide a heuris- 
tic interpretation regarding the fractional Neumann boundary conditions 
(I3.2p imposed on the parallel plates as their deviation from the perfect con- 
ductivity and infinite permeability. 
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4. Low and High Temperature Expansion and Limit of the Free 

Energy Density 

In this section, we consider the low and high temperature limits of the 
free energy density. For this purpose, a generalization of the Chowla-Selberg 
formula for Epstein zeta function (see e.g. |614 I62j ) is particularly useful. 
We have 

„ , , . cos(27rn/ii) 

Z B ( S ;ci,c 2 ;0;h)-2^ 



-J (an 2 ) 8 



(it 



2 i ~* J poo ^ 

= T(s) ^ cos(27rmh 2 ) / t 3 ' 1 £ e 

^ ' m=l n=-oo 

— OO pQQ OO 2 

- J^ F cos(27rm/ i2 ) / t s ^ V e^ 2 ™ 2 ^^ 1 ^. 

Here we have used the Poisson summation formula. If hi = 0, then we have 
to separate the n = term and obtain 

(4.1) 

/ x 2-v/tt r (s — h) cos(27rm/i9) 
Z £ ( S ; Cl , c 2 ; 0; h) = 2^(2*) + J (1/2 / 2 / £ ^ ^ 

c 2 V Cl 1 V S J m=l 

^ oo oo , x s-(l/2) / / \ 

8i/7T r—K ^— «, / ITU \ / /C2 \ 

+ \ >^ >^ cos 27rm/i 2 ^»-(l/2) 27r A/ — mn • 

If < hi < 1, then 
(4-2) 

„ . , s . cos(27rn/ii) 

Z E (s; Cl ,c 2 ;0;h) = 2q s £ ■ 



n 2s 
n=l 



+ — > > cos 27rm/i 2 — K s-(i/2) 2tt, / — m n - hi . 

4.1. Low Temperature Expansion. 

By taking c\ = l/ai,c 2 = l/a 2 ,h = (r?/2,0) in (|4TTj) and (|4~2j) , we have 
the low temperature (T < 1 or ^ < 1) expansion of the free energy density 
4323|), i.e. when rj = 0, 

= - + + + * x ,„) dg+ i r( ^i) — g J fa( g ) 

2 (D/2) + 3^+(3/2) e D/2 - - D/2 / nm \\ 

V 2 / m=l n=l s ' / 
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and when 77 7^ 0, 



oD , 2 3 (°+ 1 ) 



D+l 



E 



cos(7rnr/) 



D+l 



+ 



2 (D/2)+2 7r 2D+(3/2)^D/2 



m=l 
x x / 1 

E E 



77 

n 

2 



D/2 



D/2 



77 

n 

2 



d D+i r (D+i) 

\ 2 J m=l n=— 00 

From [B3], pg 223 , we have the following asymptotic expansion for K D / 2 ( Z ) : 



(4.3) 
where 



K u (z) 



^ (MO 



(^ fc )= 2 i!n( 4 - 2 -^- 1 ) 2 )- 



i=i 



When 1/ is equal to half of an odd integer, the sum in f|4.3[) is finite and the 
right hand side of (|4.3p is the exact formula for K u (z). From the asymptotic 
expansion (|4.3p . we see that when £ — > 0, 



EE 

m=l n=l 



D/2 



if 



D/2 



ran 

T 



is exponentially decay and the leading term is obtained by setting m = n 
1, which results in 



00 00 



EE 



D/2 



K 



D/2 



nm 

T 




2 fc 



m=l n=l v ' \ k=l 

Similarly, when < 77 < 1, the m = 1, n = term gives 



00 00 



77 

n 

2 



^ ^ (m 

m=l n=— 00 x 

V fe=i 



D/2 



77 

n 

2 



I 



CD/2, fc) 

7] k 



i k + 0[e 



mm{ 17 ,l-(r ) /2)} 

£ 



When 77 = 1, the m = 1, n = ±1 terms give 

D/2 



X X , j 

^ ^ (m 



m=l n=— 00 
1 



n 



2 (D-2)/2 



-f^D/2 

£(l + J3(/3/2,fc)£*)e-£+o(e-f). 



fc=i 
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These imply that for low temperature T< 1, when rj = 0, 
(4.4) 

, m.lU/Un 7T D / 2 r(f) D 



when < Tj < 1, 
(4.5) 

^ X , M « / £ (jgj ^ cos(7rn?7) 
Ja ' x ^ d D I 2 Z3 + 1 vr( D + 1 )/ 2 ^ n^+i 



m=l 

nn , _ * \ \ 

min{?7,l-( 



+ (W^ H |(^),^ +o| , 



and finally when 77=1, 

D+1 



(4.6) /„,,„ = - ^(- JJJ D J 1)/2 (1 - 2-") Cr (D + 1) 



e 2 « + O [ e « 



2 D 

\ k=l / / 

From these, we also find that the zero temperature energy density is 

f0 a x,» a r (z£) V cos(7rn?7) 

•>a;x^ d D 2^+i 7 r( D + 1 )/ 2 ^ n D + l ' 

m=l 

This term depends on r/. For D = 2, 3, 4, 5, the relation between the nor- 
malized zero temperature energy density d D f®. j) and rj is shown in Figure 
1. When rj = 0, its value 

f o _ <W r (^) , m ■ n 
• / a;x,M 2 D + 1 7r(' D + 1 )/ 2 1 ; 

is negative, and when rj = 1, its value 

f o _ q x ,„a rg) "(-!)» 
•>a;x,^ d D 2 d+i 7T ( d + 1 )/ 2 ^ n D+l 

n=l 

r ( D + 1 \ 

- (1 " 2 } SB 2 W# + D/2^( Z) + 1 ) 
is positive. We are going to show in the Appendix that the function 

(4.7) <B n (x) = ^ ^— -L , 0<x<l, n>2 

k=l 
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0.04 




0.2 0.4 0.6 0.8 1 

eta 



Figure 1. The normalized zero temperature free energy 
density d D f®. rj for D = 2, 3, 4, 5 and a = 1. The horizontal 
axis is the r] axis. 

is increasing in the interval [0, 1/2]. Consequently, when r\ changes from 
to 1, the zero temperature energy density increases, and it changes from 
negative to positive, so that the nature of the force in the system changes 
accordingly from attractive to repulsive. For a specific D, there is a transi- 
tion value so that the force is attractive when ij £ [0,i]£>) and the force 
is repulsive when r\ £ (t]d, 1]. We tabulate some values of r\r> in Table 6.1. 



Table 6.1 



D 


VD 


D 


VD 


2 


0.4226 


3 


0.4617 


4 


0.4807 


5 


0.4902 


6 


0.4951 


7 


0.4975 



From this table, we see that rj£> is an increasing function of D. We have 
verified numerically that this is true for all D < 45. On the other hand, we 
can in fact show mathematically that t)d < 0.5 for all D (see Appendix). In 
Figure 2, we show the graph of rjo as a function of D. 

From (|4,4p . (|4.5p and (|4.6p . we also find that when x / |Uor (x, ^) = (0, 0), 
the thermal correction to the zero temperature energy decays exponentially, 
whereas if \ = A 1 7^ 0> there is a term proportional to T D . 

When £ <C 1, the dependence of the normalized free energy density on £ 
and f] = n — x f° r D = 3 is shown in Figure 3 and Figure 4. 

4.2. High Temperature Expansion. 
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Figure 2. rj D (eta D ) function of D. 




Figure 3. The normalized free energy density d f a -,o,r) f° r 
D = 3 and a = 1. The a; and y axes are the £ (xi) and 
(eta) axes. 



Take ci = \jai,oi = l/ai,h = (0, 77/2) in (14. ip . we have the high tem- 
perature (T S> 1 or £ S> 1) expansion of the free energy density (I3.23|) . 
.i.e. 

^adT(^) ( ? (2**Z\ D+ \ (D 1 H 1 ^ D+ ^r(f) - eos(^) 

9 (D/2)+4„2D+(5/2) oo oo /rr)N D/2 \ 

+ d n»rW) ;'° +2>/2 EE ^.("V) V 

V 2 / n =lm=l " / 
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Figure 4. The normalized free energy density dP f a -,o,r) for 
D = 3 and a = 1 when £ = 0,0.05,0.1,0.15,0.2 respectively. 

Using the asymptotic expansion of the modified Bessel function (|4.3f) . we 
find that if T > 1 (or equivalently £ > 1), 

(4.8) 



(O-D+DAr f £±1A , fjD , u+d+i , r (f) f v cos(7rn??) 




n=l 



n- 1 



The leading term 

is proportional to T £)+1 , and is independent of rj. When D = 3, it gives 

f OO,l _ _ 71-2 ^ rpi 

which is called the Stefan-Boltzmann term. The next leading term of the 
energy density at high temperature is proportional to T, with proportional- 
ity constant depends on rj. The rest of the terms decay exponentially. From 
this, we can conclude that when the temperature is large enough, the effect 
of different boundary conditions is not significant and the system exhibits a 
universal behavior regardless of the boundary conditions. When £ 3> 1, the 
dependence of the normalized free energy density on £ and tj = fi — x f° r 
D = 3 is shown in Figure 5. 

As we explain in Section [3T6l if we take ct XjA1 = 2, uj x ^ = 0, a = 1, rj = 
and D = 3 in the energy density fa\x,^ (|3-23j) . we obtain the Casimir energy 
for electromagnetic fields confined between perfectly conducting parallel in- 
finite plates (|3.24p . Therefore, we obtain from (J44]) and (|4.8p the low and 
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x 10 5 




Figure 5. The normalized free energy density d f a -,o,r) f° r 
D = 3 and a = 1. The x and y axes are the r] (eta) and £ 
(xi) axes. 



high temperature limit of the Casimir energy (j3.24j) : 

fca, + yCk(3){ 3 + » 2 ft 2 + C) e-i ) + O (e-i) , f « 1, 

fe. ~ - ^ + |fe(3) + (^f + 1) + o ( e - 8 - 2 «) , « » 1, 

agree with the result of [60J. On the other hand, if we take o x ^ = 2, 
^x.m = 0, a = 1, 7/ = 1 and D = 3, we obtain the Casimir energy for 
electromagnetic fields confined between one perfectly conducting and one 
infinitely permeable parallel infinite plates. Therefore, from (|4.6p and (|4.8|) . 
we find that the low and high temperature limits of the Casimir energy 
density of this system are 




£< i, 
£» 1. 



These agree with the results in |33| . 

5. Temperature Inversion Symmetry 

Since the observation of the symmetry between low and high tempera- 
ture exhibited by the Casimir energy between perfectly conduction parallel 
plates (|3.24p pointed out by Brown and Maclay in [60], there have been a 
number of papers devoted to the discussion of the temperature inversion 
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symmetry of different systems [331 IM1 [35], [361 EH EH1 US]. Here we want 
to point out the mathematical origin of this symmetry, and show that in 
some particular cases, the free energy density (|3.22p exhibits temperature 
inversion symmetry. 

Observe that when h = = g, the Epstein zeta function (|3.10p is com- 
pletely symmetric with respect to a\ and a 2 - In particular, if we define 

H s (w)= V ' — s = Z E {s;l,w 2 ;0,0), 

, , ^ nf + (n 2 w) 2 ) 

{n 1 ,n 2 )ei, 2 v 1 v ' 1 

then 

Z E (s; ai ,a 2 ;0,0)=a- s H s ( *[^) = a^ s H s ( x [^ 

VV a iJ VV °2, 

The symmetry of the Epstein zeta function expressed using H s is the relation 

(5.1) H s (w)=w- 2s H s (- 

\w 

In our case, \fa^Ja\ = 2tt^ = 2dT and formulas in the form (|5.ip precisely 
gives a relation between low and high temperature. 

Using the formula (|3.23[) for free energy density, when x = /i, we have 
g = and therefore 

a U;x,x ~ „ „ 3(p+i) a 2 H s±i [ \ — ) ~ a x,n u x,v 9 * CM^J- 

2 D+2 7T 2 2 VV a l/ 1 

The second term is zero except when x = ^ = or 1. The first term, 
denoted by JFg , is a function of £ and is equal to 

Z 2 

From (|5.ip . it satisfies the inversion symmetry 

1 



4vr 2 £ 



For D = 3, a = 1, this is precisely the symmetry observed in [601 [36 J for 
electromagnetic field confined between parallel perfectly conducting plates. 
Therefore, when x = ^ / 0,1, the normalized free energy density of a 
massless fractional Klein-Gordon field confined between two parallel hyper- 
planes dP fa;x,n nas a complete temperature inversion symmetry. When 
X = [i = 0, 1, the symmetry is broken by a term proportional to £ d Cr(D). 

When 7] = ±1 or equivalently, (x, f^) = (0,1) or (1,0), from (13.23H we 
have 

d D fa;0,l = ~ - 3( D+ \) (V a l a 2GQ D+1 Y't 2~T Wd+IW 
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We can rewrite the double sum as 



V v , (~l) n 
^ ( ai m 2 + a 2 n 2 )( D + 1 )/ 2 

m£Z neZ 



Therefore, the normalized free energy density (i D / a; o,i can be written as a 
sum of two functions in £, iFi(^) and 2^1 (C) where 



^(0=fr(^±i)vr^ +1 ^(2<). 



Using (|5.ip . we find that each of these functions satisfies an inversion sym- 
metry 



When Z? = 3, a = 1, this is what observed in [33] for electromagnetic field 
confined between parallel plates under Boyer's setup. 

For generic 77, there was no temperature inversion symmetry since the 
components in g are not symmetric. However, for some particular rational 
values of 77, we can use the same trick as in the case 77 = 1 and write the 
normalized energy density as a sum of a few functions such that each of 
them has temperature inversion symmetry. For example, when 77 = 2/3, 
using the fact that when f(x) is an even function, 



cos ( — ) f(n) 

mod 3 



^e™"/(n)=/(0)+2 Yl /(»)+ E 

riGZ \n>l,n=0 mod 3 ri>\,n=\ i 

+ E cos (f) /(»)) 

n>l,n=2 mod 3 ^ ' / 
oo 

=/(o)+3 /h-E/w 

n>l,n=0mod3 n=l 



2 v ' 2 
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we can write d f a - x x ±(2/3) as a sum of two functions ^F^iOi 2^*2(0 given 
by 

iMO = - (^±1) ^e +1 Ho^ (6*0 , 

^(0 =fr(^±i)vr^^ +1 ^(2vre). 
Each of these functions satisfies temperature inversion symmetry 

im) = (6^) D+1 (^) , «o = (2<) D+1 & (^) . 

6. Conclusion 

We have introduced a new type of boundary condition called fractional 
Neumann condition which involves vanishing fractional derivative of the field 
in the study of Casimir effect of fractional massless Klein-Gordon field con- 
fined between a pair of parallel plates. By imposing this fractional Neumann 
conditions on the plates allows the interpolation between the usual Dirichlet 
and Neumann conditions. Our results indicate that there exists a transition 
value for the difference between the orders of the fractional Neumann condi- 
tions in the two plates for which the Casimir force changes from attractive 
to repulsive (or vice versa). It is interesting to note that for sufficiently 
high temperature, the Hemholtz free energy density is dominated by a term 
independent of boundary conditions. Conditions for temperature inversion 
symmetry to hold are also discussed. 

We would also like to point out that despite a few decades of work on tem- 
perature dependence of Casimir effect, there still exist debates on this topic. 
The main issue of the recent controversy lies in the thermodynamic con- 
sistency of the computed Casimir force between real metals and the Drude 
dispersion relation (see references [631 ESI ESI E2 EH] for both sides of the 
controversy). It has to do with the controversy of inclusion/exclusion of 
the TE (transverse electric) zero mode. Some authors [691 EE \U\ claimed 
that the Drude relation does not provide a consistent explanation of re- 
cent experimental results, in particular it is in conflict with the Nernst 
theorem. They proposed to replace the Drude relation by the plasma re- 
lation. On the other hand, Hoye, Brevik, Aarseth, Ellingsen and Milton 
[651EZ1IZ21E31E31E51E61E71E8] have argued in favor of the exclusion of the 
TE zero mode. They have derived analytical results using Euler-Maclaurin 
formula, which in the limit T — * are consistent with the Nernst theorem 
|79j . They have also carried out numerical calculation of the free energy 
and obtained results which agree with analytic results to a high degree of 
accuracy. These authors also proposed an experimental setup to test such 
results |78j . We plan to discuss in detail the Casimir energy of fractional 
electromagnetic field and the issue of inclusion/exclusion of the TE zero 
mode in a future work. 
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Finally, we would like to suggest some other possible directions for further 
work. The extension of our discussion to a p-dimensional cavity embedded 
in D-dimensional space, with p < D is currently under consideration. How- 
ever for the generalization of the above results to non-flat space is likely to 
encounter highly non-trivial mathematical problems since one needs to deal 
with fractional operators in curved space. Another interesting generalization 
involves fractional Klein-Gordon field with fractional Neumann boundary 
conditions of variable fractional order, which allows variable Casimir energy 
or force at different point in space. Such a problem again requires results 
from derivatives and integrations of fractional variable order, a subject which 
is still at its infancy. 
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Appendix A. The Function <B n (x) 
In this appendix, we are going to show that the function (|4.7p 

cos(2-7rA;x) 



k=l 



k h 



is increasing and has exactly one zero in the interval [0, 1/2]. We are also 
going to show that this unique zero is less than 1/4. 

First, we show that ^Sn^x) is increasing and has exactly one zero in the 
interval [0, 1/2]. As a matter of fact, for n even, the function 2$ n (x) is well 
known. From 9.622 of [57], we have 

(27r) 2n 

(A.l) <B 2n (x) = (-ly-iLJL.Btoix), 0<x<l, 

where Bk(x) is the fc-th Bernoulli polynomial defined by 

te tx 



= B k {x) k 

e* - 1 ^ k\ 

k=0 

The explicit formula for Bk{x) for 1 < k < 5 is given in Table A.l. 



Table A.l 


k 


B k {x) 


1 


r> X 2 , 1 


2 


x z - x + I 


3 


/y>3 3 '2 _i_ _1 , r 

•X- 2 I 2 


4 


-r 4 — 9-1-3 _i_ J2 1_ 

*ju ^ \ uij 30 

rJO 5.0^4 1 5™3 ±.rp 

lAJ 2 1 *^ «A> ^ iXj 


5 
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It is well known that for all k > 1, B' k+l (x) = (k + l)B k {x). From 9.622 
of [57] again, we have 

(A.2) 

_ (-l) n 2(2n - 1)! ^ sin(2fc7rx) (k = 1, 0<x<l 

Jhn-lW- (27r) 2„-l A, fc 2n-l ' jjfc>2, < X < 1 ' 

From this it is easy to verify that B 2n (x) is increasing and has exactly one 
zero in the interval [0,1/2] (see e.g. [80J). For convenience, we repeat the 
argument here. As is easily verify from (jA.lj) . 



ya-l (gn)l 

(2tt)< 



B 2n (0) =2(-l)"- 1 A-^CK(2n) 



(2n)\ 



B 2n (l/2) =2(1 - 2 1 - 2 «)(-l)«^ r Ci ? (2n), 

which shows that -B 2 n(0) and -B 2 n(l/2) has opposite sign and is nonzero. 
It also implies that B2 n {x) must has at least one zero in [0,1/2]. On the 
other hand, we find from ([A~2|) that B 2k+1 (0) = B 2k+1 (l/2) = for all 
k > 1. Now if for some j > 1, B 2 j(x) has two zeros in [0,1/2], then its 
derivative 2jB 2 j_\{x) has a zero in (0, 1/2). Since -E>2j-i(0) = B 2 j-i(l/2) = 
0, this in turn implies that its derivative (2j — l)B 2 j~ 2 (x) has two zeros in 
[0, 1/2]. Continuing this argument, we find that B\(x) must have a zero in 
(0, 1/2). This gives a contradiction since B\(x) = x — (1/2) does not have 
any zero in (0,1/2). This shows that B 2n {x) has exactly one zero in the 
interval [0, 1/2] and B 2n _\(x) does not have any zero in the open interval 
(0, 1/2). The latter implies that B 2n _i(x) must be either always nonnegative 
or always nonpositive in the interval (0, 1/2). Therefore, B 2n (x) is monotone 
in [0, 1/2]. This completes our argument for 53 2n (x). 

To verify the statement for 2? 2n __i(x), n > 1, we define the functions 

~ / \ sin(2A;7rx) f n = 1, < x < 1 

2>n(x) = E 



k n \n > 2, < x < 1 

k=l \ 



and let 



(-l)™- 1 

C *n-l{x) = (2 ^ 2w _i %2n-l(x), C7 2 „(:r) = ( ^ )2n £>2n(a) 

for all n > 1. Then it is easy to verify that C^ +1 (x) = C n (x). Moreover, for 
all n > 1, 



(2tt)= 



^n + i(0) fcrfan + l), 



f— 1 "i n 

C 2n+1 (l/2) = - (1 - 2- 2 ")^-pA_C«(2n + 1), 
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C 2n (0) = C 2 „(l/2) = 0. On the other hand, for < x < 1, 

oo 2,irikx / 

53i(x) + iDi(x) = ^ = " log (1 - e 2nix ) = - log(2 sin(vrx)) - ivr f x - 

k=i ^ 

Therefore, 

d(x) = ^-53i(x) = -— log(2sin(7nr)) 

and it is easy to verify that C\{x) is decreasing on (0,1/2), positive on 
(0,1/6), negative on (1/6,1/2] and zero at x = 1/6. Since C' 2 {x) = C±(x), 
we find that C 2 is strictly increasing on [0, 1/6] and strictly decreasing on 
[1/6, 1/2]. Since C 2 (0) = C* 2 (l/2) = 0, C 2 {x) > for all x £ (0, 1/2). The 
same argument used for £? 2n then shows that C 2n+ i(x) is monotone and 
has exactly one zero in the interval [0, 1/2], thus verifying the statement for 

»2n-l(x). 

Now since <B n (0) = Oi(n) > 0, Q3 n (l/2) = -(1 - 2 1 - n )( R (n) < 0, to show 
that the unique zero of 25 n (x) is less than 1/4, it is enough to show that 
25 n (l/4) < 0. A straightforward computation gives 

/l\ 1111 1 /1\ 

<Bn \ v 4y~"F + 4^"6^ + 8^ + '"~ 2» n \2J < ' 

verifying our claim. 
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